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' Path integral techniques and Green functions formalism are applied to study the (time) temper- 

I ature dependent scattering of a polaronic quasiparticle by a local anharmonic potential in a bath of 

diatomic molecules. The electrical resistivity has been computed in any molecular lattice dimension- 
' ality for different values of electron-phonon coupling and intermolecular forces. A broad resistivity 

peak with non metallic behavior at temperatures larger than ~ WOK is predicted by the model for 
sufficiently strong polaron-local potential coupling strengths. This peculiar behavior, ascribed to 
purely structural effects, is favoured in low dimensionality. 
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4L ! I.INTRODUCTION 

. A considerable amount of theoretical work has been 
devoted to investigate the condition of polaron forma- 
tion and the polaronic features in real materials 
V ' While there is growing evidence that fundamen- 
I tal properties such as the polaron size, effective mass and 
^ ' ground state energy are essentially similar in any dimen- 
O ' sion [jll],^, it is still unclear to which extent transport 
properties in polaronic systems (markedly the electrical 
resistivity behavior versus temperature) depend on the 
CN , lattice structure and dimensionality. Besides being con- 
^ ' ceptually relevant this question has become actual in 
connection with the discovery of unusual effects in un- 
derdoped high Tc superconductors. Infact the presence 
of local lattice distortions with polaron formation has 
been envisaged in the high Tc systems and signs of en- 

■ hanced anharmonicity for some in- and out of plane oxy- 
, gen modes have been detected in underdoped compounds 

■ by several groups In this paper we focus on the 
problem of the interaction between a polaronic quasipar- 

' tide moving through a molecular lattice and a local struc- 

' . I tural instability modelled by a double well potential in 

^ . its two state configuration. In particular we derive the 

Q effective interaction strenghts arising from this peculiar 

Q . scattering mechanism and we study the effects on the 

ILjl ■ electrical resistivity both of the intermolecular forces and 
of the lattice dimensionality. 
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II. THE MODEL 



Our analysis starts from the following r dependent 
Hamiltonian where r is the time which scales as an in- 
verse temperature according to the Matsubara Green's 
function formalism: 

Hoir) - e(.g)£t(T)a(r) + ^ Wqa]j(r)aq(r) + HrLsir) 

q 



1 



Q{t)^ ~Qo + —{t-U). (0.1) 

To 

Ho{t) is the free Hamiltonian made of: a) a polaron 
created (distroyed) by c' {t){c{t)) in an energy band e{g) 
whose width decreases exponentially by increasing the 
strength of the overall electron-phonon coupling constant 
g, e{g) = Dexp{—g'^); b) a lattice of diatomic molecules 
whose phonon frequencies Wq are derived analytically for 
a linear chain, a square lattice and a simple cubic lattice 
through a force constant model; c) a local anharmonic 
potential shaped by a Two Level System (TLS) in its 
symmetric ground state configuration. Q{t) is the one di- 
mensional space-time hopping path followed by the atom 
which moves between two equilibrium positions located 
at ±Qo- To is the bare hopping time between the two 
minima of the TLS and U is the instant at which the ith- 
hop takes place. One atomic path is characterized by the 
number 2n of hops, by the set ofU {0 < i < 2n) and by Tq. 
In the last of eqs.(l), we assume that the class of r-linear 
paths yields the main contribution to the full partition 
function of the interacting system. The closure condition 
on the path is given by: (2n — 1)ts + 2nro = (3, where /3 
is the inverse temperature and Ts is the time one atom is 
sitting in a well. The interaction is described by Hintir) 
with A being the coupling strength between TLS and po- 
laron, XQ(t) is the renormalized (versus time) tunneling 
energy which allows one to introduce the r dependence 
in the interacting Hamiltonian JTst . 

Following a method previously developed llfl] in the 
study of the Kondo problem, we multiply AQ(r) by a 
fictitious coupling constant s (0 < s < 1) and, by dif- 
ferentiating with respect to s, one derives the one path 
contribution to the partition function of the system 

^^f Z{n,U) \ ^ r f\rQ{T) lim G(t,t'), 

\ ^0 J Jo Jo t'^t+ 

(0.2) 

where, Zq is the partition function related to Ho and 
G(r, T )s is the full propagator for polarons satisfying a 
Dyson's equation: 

G{t, t)s= G"(r, r ) + s / dyG°(r, y)XQ{y)G{y, r ) 

(0.3) 







The polaronic free propagator G*^ can be derived ex- 
actly in the model displayed by eqs.(l) |0. We get the 
full partition function of the system by integrating over 
the times ti and summing over all possible even number 
of hops: 
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Zt ^ ^oV" / \ —exp -PE{n,ti,To) 

tr'Jo ro Jo TO L 

PE{n,U,To)=L-{K^ + K^)f2(*-l-^y (0.4) 

with E{n,ti,To) being the one path atomic energy. L, 
which is a function of the input parameters, is not essen- 
tial here while the second addendum in eq. (4) is not local 
in time as a result of the retarded polaronic interactions 
between successive atomic hops in the double well po- 
tential, and are the one path coupling strengths 
containing the physics of the interacting system. 
(negative) describes the polaron-polaron attraction me- 
diated by the local instability and (positive) is re- 
lated to to the repulsive scattering of the polaron by 
the TLS. Computation of E(n,ti,To) and its derivative 
with respect to tq shows that the largest contribution to 
the partition function is given by the atomic path with 
Tg — Q. The atom moving back and forth in the dou- 
ble well minimizes its energy if it takes the path with 
the highest tq value allowed by the boundary condition, 
that is with {To)max — {2nKBT)^^. This result, which 
is general, provides a criterion to determine the set of 
dominant paths for the atom at any temperature. Then, 
the effective interaction strengths < > and < > 
can be obtained as a function of T by summing over n 
the dominant paths contributions: 

<K^ >= -(AQo)'s'exp(2^Aq)^^qC^2^~ 

q q 

< -f3{XQo)^B^exp{3 ^ A^) ^ A^ujH (0.5) 

q q 

with: B = {nri^) — ^)Gxp[—g'}2,^ctgh{Pujq/2)) and 

Nq is the phonon occupation factor and npi^ig)) is the 

Fermi distribution for polarons. / = 'Y^^=i{'^^)max ^-i^d 
N is the cutoff on the number of hops in a path. The 
particular form of (ro)maa; suggests that many hops paths 
are the relevant excitations at low temperatures whereas 
paths with a low number of hops provide the largest con- 
tribution to the partition function at high temperatures. 
Since the effective couplings which determine the resis- 
tivity depend on [tq)'^^^^, hence on (through /), a 
relatively small cutoff {N ~ 4) ensures numerical conver- 
gence of eqs.(5) in the whole temperature range. On the 
other hand the non retarded term L in eq.(4) has a slower 
behavior, therefore a larger cutoff should be taken 
at low temperatures where computation of equilibrium 
properties such as specific heat is strongly influenced by 
many hops atomic paths between the minima of the dou- 
ble well excitations. 

The lattice Hamiltonian is made of diatomic sites 
whose intramolecular vibrations can favor trapping of 
the charge carriers p8[ . The mtramolecular frequency 
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largely determines the size of the lattice distortion associ- 
ated with polaron formation while the dispersive fea- 
tures of the phonon spectrum controlled by the first (wi) 
and second (0^2) neighbors mtermolecular couplings are 
essential to compute the polaron properties both in the 
ground state and at finite temperatures |^o|] . The range 
of the intermolecular forces is extended to the second 
neighbors shell since these couplings remove the phonon 
modes degeneracy (with respect to dimensionality) at the 
corners of the Brillouin zone thus permitting to estimate 
with accuracy the relevant contributions of high sym- 
metry points to the momentum space summations 
Then the characteristic frequency Co, which we choose as 
the zone center frequency, is: = 10^ + zujI 4- z„„„[jj|, z 
is the coordination number and 2„„„ is the next nearest 
neighbors number. Hereafter we take frequencies values 
which are appropriate to systems with rather sizeable 
phonon spectra. 

We turn now to compute the electrical resistivity due 
to the polaronic charge carrier scattering by the impurity 
potential with internal degree of freedom provided by the 
TLS. Assuming s wave scattering, one gets [ p2| , p3[ 

p = posin^rj 

po = o— (0.6) 

7Te^vl{No/V) h 

where, Ug is the density of TLS which act as scatterers, 
vp is the Fermi velocity, V is the cell volume, Nq is the 
electron density of states, e is the charge and h is the 
Planck constant. The phase shift ij to the electronic wave 
function at the Fermi surface is related to the effective 
interaction strengths < > and < >. 

The input parameters of the model are six that is, the 
three molecular force constants, g, D and the bare energy 
AQo- Qo can be chosen as ~ 0.05A consistently with re- 
ported values in the literature on TLS's which are known 
to exist in glassy systems ^4-27 , amorphous metals 



A15 compounds and likely in some cuprate superconduc- 
tors p9|-pl||. In these systems the origin of the TLS's 
is not magnetic. The bare electronic band D is fixed at 
O.leV. 

Let's start by studying (Fig.l) the resistivity behavior 
in a 3D lattice as a function of the overall electron-phonon 
coupling g. X is fixed at 70QmeVA~^ which means a 
bare tunneling energy of ~ 35meV and the molecular 
force constants are loq — 60meV, loi — ^QmeV and 
UJ2 = 20meV, respectively. At low temperatures p(T) 
approaches the unitary limit for any of the four g values 
thus displaying the peculiar effect of the anharmonic po- 
tential with internal degree of freedom. At T ~ 150K, 
p develops a g-dependent broad peak which softens and 
finally disappears by increasing g in the strong coupling 
regime. We note that at g > 1 the polaron energy be- 
comes smaller than the atomic tunneling energy hence 
polaron scattering by the TLS is essentially diagonal. On 
the contrary, in the intermediate coupling regime (g in 
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the range 0.5 - 1) the incoming polaron can release a suf- 
ficiently high amount of energy, off diagonal scattering 
by the impurity potential prevails and a broad resonance 
peak emerges. At weaker g the polaronic picture would 
lose any validity. 

Next we turn to consider the effect of the lattice dy- 
namics on the p versus T behavior. In Figs. 2 and 3, the 
normalized resistivities are reported on for the ID and 
2D molecular lattices, respectively. Let's set 5=1 which 
ensures both polaron formation and presence of the reso- 
nance peak also in ID and 2D systems while the tunneling 
energy (as in Fig.l) is in the range of the values estimated 
by EXAFS investigations in high Tc systems with double 
site distributions for the apical oxygen atoms (2^,^ . By 
taking ujq — QOmeV and lui ~ 50meV we choose for the 
first neighbors intermolecular coupling model a charac- 
teristic phonon energy Co ~ O.leV. This choice allows us: 
i) to treat correctly the ground state polaron properties 
versus dimensionality ||21| , ii) to reproduce the sizable 
phonon energy of c-axis polarized mode due to apex oxy- 
gen vibrations coupled to the holes in the Cu-0 planes 
of Y Ba2Cuj,0-!-s- In underdoped compounds (whose c- 
axis resistivity displays the unusual non metallic behav- 
ior) this mode appears enhanced in energy . Although 
our simple cubic lattice does not account for the details 
of the structural effects of YBCO we are therefore in 
the appropriate range of parameters to capture the main 
features of the lattice polarons scattered by local instabil- 
ities in those compounds. The role of the second neigh- 
bors couplings is emphasized in Figs. 2 and 3 by varying 
UJ2 which however should not exceed toi. Increasing W2 
the polaron spreads in real space and becomes lighter. 
Accordingly the height of the peak decreases. We also 
note that the strength of the intermolecular forces infiu- 
ences the position of the peak versus temperature which 
results from a balance between competing attractive and 
repulsive interactions (eqs.(5)). At fixed parameters, the 
height of the ID peak is roughly twice as much that of the 
2D one. Comparing the dotted 2D curve {uj2 = 20meV) 
with the corresponding case {g = 1 plot) in Fig.l, one 
sees that the 3D peak is further reduced and its absolute 
value lies below the residual resistivity value. 

The dynamics of the TLS and its coupling to the charge 
carrier strongly affect the transport in any dimensional- 
ity as Figs. 4 make evident. When A is small (i.e. below 
300TneVA~^ in ID) polaron and TLS are weakly cou- 
pled, off diagonal scattering is unlikely to occur and the 
conductivity is metallic-like. Increasing A, the atomic 
tunneling energy becomes of order of the polaron energy 
and the conditions for resonant scattering are established. 
Note that by varying A the resistivity peak does not shift 
versus temperature while the A threshold value for the 
appearance of metallic conductivity is larger in higher 
dimensionality. This means that 3D systems can sustain 
a larger degree of anharmonicity (than low dimensional 
systems) and still exhibit metallic transport properties 
while trapping of the charge carrier by the anharmonic 
impurity potential is favoured in ID. It is this trapping 
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which causes the broad resistivity peak and non metallic 
transport at T larger than ~ lOOif in ID. 

We investigate further the effect of the molecular forces 
(Figs. 5) by increasing the miramolecular energy and 
decreasing the miermolecular energies with respect to 
Figs. 4. The ID peak (compare Fig. 5(a) and Fig. 4(a)) 
is reduced for any A value whereas the 2D and 3D peaks 
are here larger than in Figs. 4(b) and 4(c) respectively. 
The reason is the following: in ID, the intramolecular 
coupling is determinant because of the low coordination 
number; as a result, by enhancing Wq the carriers be- 
come lighter and the absolute resistivity is reduced. On 
the contrary in 2D and mostly in 3D, long range effects 
are more effective hence, the stronger wq (with respect 
to Figs. 4) is more than balanced by the weaker and 
UJ2', as a result, the characteristic 2D and 3D phonon fre- 
quencies are smaller, polarons become less mobile and 
the resistivity peaks are accordingly higher. 



III. CONCLUSION 

The path integral formalism has been applied to study 
the time retarded interacting problem of a polaron scat- 
tered by a local anharmonic potential in a lattice of 
diatomic molecules. I have derived the full partition 
function of the system and obtained analytically the ef- 
fective coupling strengths as a function of temperature. 
The electrical resistivity has been computed in any lat- 
tice dimensionality for a large choice of input parame- 
ters. Strength of the overall electron-plionon coupling 
{g), strength of the polaron- local potential coupling (A) 
and strength of the molecular forces interfere, giving rise 
to a rich variety of resistivity behaviors versus tempera- 
ture. As a main feature, when the conditions for resonant 
scattering between polaron and local double well poten- 
tial are fulfilled, a broad resistivity peak shows up in 
the 100-150K range. However, the shape of the peak es- 
sentially depends on the dynamics of the local potential 
(tuned by A) and, for sufficiently low atomic tunneling 
energies, metallic conductivity conditions are recovered. 
Generally, in 3D the resistivity maximum is less pro- 
nounced than in low dimensionality and, for sufficiently 
strong intermolecular couplings, the height of the peak 
lies below the residual resistivity value, being therefore 
hardly visible in experiments. On the other hand, in 
ID and 2D, the maximum can easily become larger than 
the residual resistivity by a factor ~ 3 — 4 for realistic 
choices of tunneling energies and lattice force constants. 
Then, non metallic resistivities can be expected in low 
dimensional polaronic systems with local lattice insta- 
bilities at least for temperatures above ~ lOO-ftT. As a 
final remark one should add that, while the intra- and 
mfermolecular energies have been varied throughout the 
paper as apparently independent parameters, their val- 
ues for real systems can be obtained (for instance) by a 
least squares fitting procedure to experimentally known 
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physical quantities Q. 



FIG. 1. Electrical Resistivity normalized to the residual 
(T — 0) resistivity for four values of electron-phonon cou- 
pling g. The bare TLS energy XQo is 35meV. The force con- 
stants which control the phonon spectrum in the 3D lattice 
are: cjq ~ 60meV, u)\ — 50me V, 0^2 = 20meV 



FIG. 2. ID Electrical Resistivity as a function of the second 
neighbors intermolecular force constant, g — 1. XQq, ujq and 
uji are as in Fig.l 
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FIG. 3. 2D Electrical Resistivity as a function of the second 
neighbors intermolecular force constant. 5=1. AQo, and 
wi are as in Fig.l 



FIG. 4. Electrical Resistivities for six values of the 

poIaron-TLS coupling A in units meVAr^. g = 1- 
Wo = 60meV,uJi = 50ineV,uJ2 = SOmeV. (a) one dimen- 
sion; (b) two dimensions; (c) three dimensions 
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FIG. 5. Electrical Resistivities for sbc values of the 
polaron-TLS coupling A in units meVA~^. g = 1- 
ujo = 80meV,uJi = 4QmeV,uJ2 = 20meV. (a) one dimen- 
sion; (b) two dimensions; (c) three dimensions 
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